Encoding relativistic potential dynamics into free evolution 
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We propose a method to simulate a Dirac or Majorana equation evolving under particular poten- 
tials with the use of the corresponding free evolution, while the potential dynamics is encoded in a 
static transformation upon the initial state. We extend our results to interacting two-body systems. 
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I. INTRODUCTION 

The last years have witnessed an increasing interest in 
simulating dynamics coming from the relativistic quan- 
tum mechanics realm in different physical systems, e.g. 
in trapped ions IJ, optical lattices 2, 3J, and quantum 
photonics [4,. Striking theoretical predictions related 
with the Dirac equation [S] like Zitterbewegung and Klein 
paradox [6] have been observed in these simulations. In 
particular, the proposal of simulation in one trapped ion 
[7j of the free single-particle one-dimensional (ID) Dirac 
equation has been successfully implemented in the lab 
[8]. This is also the case for the single-particle ID Dirac 
equation with some external potentials under which the 
particle exhibits Klein tunneling [5], where the experi- 
mental implementation involved two ions 110]. 

In reference [TT], it was proposed that the free Ma- 
jorana equation [12l [13] and unphysical operations like 
complex conjugation, charge conjugation, and time re- 
versal can also be simulated with two trapped ions. Be- 
sides, two-body Dirac equations have been the subject 
of recent theoretical research [lllll]. In general, simula- 
tions of single-particle equations with external potentials 
or many-body equations for interacting systems seem to 
be much more demanding than free equations. 

In this work, we show that the same setups employed 
for the simulations of the free single-particle Dirac and 
Majorana equations can also be used for simulations of 
these equations with the addition of a broad class of po- 
tentials. This is based in the following concept, which 
is the main result of this paper: any state which is a 
solution of the Dirac or Majorana equation with one of 
these potentials can be related through a static transfor- 
mation with a solution of the free corresponding equa- 
tion. Accordingly, in order to simulate the dynamics of a 
given state under certain potential, which could be cum- 
bersome in some situations, all that is needed is to ini- 
tialize the appropriate state and let it evolve under the 
free equation, which is often easier to implement. The 
method includes scalar and spinorial position-dependent 
potentials, different in the Dirac and Majorana equa- 
tions. In general, the transformation does not leave the 
probability density unchanged but this happens in some 



particular cases, showing us an additional interesting fea- 
ture, namely, under certain potentials the particle be- 
haves as a free particle. In other cases, we show that the 
method works approximately in some regions of space 
or in some parameter ranges. In particular, a massive- 
particle dynamics can be simulated in this way with a 
massless-particle equation. Finally, we show how to ex- 
tend our results to two-body situations. 

Although our method works for any dimension and 
representation, we focus here in the ID case which has a 
direct experimental connection with trapped-ion experi- 
ments simulating relativistic quantum dynamics. We use 
also the particular representations employed in these ex- 
periments. The techniques developed in this work are 
also valid for 3D but the particular results obtained in 
that case must be the focus of further research. 



II. ONE PARTICLE SYSTEMS 

We will consider the following ID Dirac equation in 
natural units {h = c = 1), 



iip = -ia^ip' + [azm + V{x)]ip, 



(1) 



where ' and ' denote time and space partial derivatives 
respectively. In the same representation, the Majorana 
equation reads 

iip = —iox ip' — i<Jy m^p* + Vix)ip. (2) 

Moreover, in ID a general potential can be written as [9] 

V{x) = fii^) + f2{x)cr,, + f3{x)<Jy + f4{x)ax. (3) 

Note that in the case of Eq. (fTl) , we could consider mass- 
like potentials taking f2ix) = m. 

A. Majorana equation 

We start from Eq. ^ where the potential is given by 
Eq. pi. We will analyze a set of sufficient conditions 
under which a transformation of the form 



ip — u{x)4>, 



(4) 



where 
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U{x) 



-i F\{x)Gx—i F2{x)ay—iF-i{x)(7z—iF4^{x) 



(5) 



can convert Eq. (pi) into the corresponding free Majo- 
rana equation for (j). (Notice that e'^ e^ ^ e^~^^ unless 
[A, B] = 0, which in general is not the case here.) To this 
end, we first notice that after applying Eq. Q the LHS 
of Eq. ([2| becomes 

iU{x)4', (6) 

while the RHS transforms into 

-io-x [U{x)(j)'+U{xy(l3] -imayU{x)*(l)*+V{x)U{x)<j), (7) 

We now impose 

iaJA{x)' = V{x)U{x). (8) 

In the cases where the derivative of the exponent of U 
commutes with U, i.e., restricting U in Eq. ((sl) to 



U{x) 



-iFj {x)aj—iF4(x) 



(9) 



with CTj = crxjCTyjaz for j = 1,2,3 respectively, we can 
write Eq. ^ as 

a, [F;(x)fT, + F:^{x)]Uix) = V{x)U{x). (10) 

This amounts to removing the potential from the dy- 
namical equation for choosing properly the relation- 
ships between the functions Fj{x) and fj{x) in Eqs. (IS]), 
(Is]). Specially interesting is the case 

F[{x) = Mx), Fi{x)^U{x), 

F2{x) ^ F^ix) ^ f2{x) ^ Mx) ^ 0, (11) 

that represents an electromagnetic potential acting on 
the Majorana particle. Under this condition, we will now 
set the commutation rules of U with the Pauli matrices. 



axU = Uax, and ayU* = U<jy, 



(12) 



while the first is simply the definition of lA, the second 
gives U* in terms oiU, i.e. 

U = aJJax, and U* = ayUuy, (13) 

which implies that Fi has to be real and F^ imaginary. 
In this case, the Majorana equation (l2| becomes 



iUcj) = lA{—iax4>' — imay(j)*), 



(14) 



where hi =IA and (j) satisfies a free Majorana equation. 

That is, if we want to simulate the dynamics of a par- 
ticle in the state ipix, t) under the ID Majorana equation 
with a potential of the form 



V{x)^hix) + U{x)a,, 



(15) 



we only have to prepare the initial state 0(x, i = 0) which 
is related with tp{x,t = 0) through 

^(x,i = 0) = e-*-^i(^)'"==-^-^''(^) (/)(x,t = 0), (16) 



and then let the system evolve under the free Hamilto- 
nian of Eq. (fl4 1 . 



Since Fi is real and ^4 imaginary, we have the following 
relationship between the probability densities: 



mx,t)\'- 



^ g-2jF4(x) 



|<^(a;,t)|' 



(17) 



Accordingly, the probability density observed for can 
be easily related with the simulated one for ip. Analogous 
relations can also be derived, for instance, for expectation 
values of observables. 

Of particular interest is the case in which F4{x) = 0, 
where the potential does not include the term with f4{x), 
because the probability density is the same for ip and (p. 
Therefore, the probability density of a state ^'(2;, t) under 
the potential V{x) = fi{x), is always the same as the one 
of the transformed free state (j){x,t). We shall illustrate 
this with an example below. 



B. Dirac equation 

We now analyze the case of the ID Dirac equation of 
Eq. ([I]), assuming that we have a massive system m 7^ 0. 
Using the same techniques as in the Majorana case, we 
can eliminate a potential of the form 



from Eq. ([l|, with the transformation 

where 

F'^ix) = h{x). 
The resulting dynamical equation for (j) is 

icf) = —iax<t> + <^z'ni4>, 
and the relation between 1-01 and |0| is given by 



IV'I 



^21mlFi(x)] 



101, 



(18) 



(19) 



(20) 



(21) 



(22) 



with Im[i^4(x)] the imaginary part of F4{x). Note that 
with this method we are allowed to simulate the pres- 
ence of complex potentials in the Dirac dynamics through 
the Hermitian free equation (21 1. This may be useful to 



study weak dissipative processes. 



C. Non-relativistic limit 

The case of a Dirac potential with only /4 7^ 0, as 
is the case for the minimal coupling with an electro- 
magnetic field, can be solved exactly and is related to 
the U(l) gauge invariance of Quantum Electrodynam- 
ics (QED). This also shows the partial relationship be- 
tween our transformations and the standard gauge trans- 
formations of quantum field theories. This kind of po- 
tential gives rise, in the non-relativistic limit p <^ m, to 



a ID Pauli-Schrodinger equation (without scalar poten- 
tial), with Hamiltonian H = [p — fi{x)Y' /2m. Since the 
non-relativistic limit amounts to the dispersive limit of 
the corresponding simulating trapped- ion equation [7], 
this shows that our method can also be applied to the 
simulation of non-relativistic dynamics. In the QED case, 
the vector potential considered in this example would be 
irrotational, such that it does not contribute to the in- 
teraction. However, we consider here a more general for- 
malism, in which an arbitrary potential V{x) = fi{x) Ux 
can be eliminated by the phase transformation. 

As we have seen in Eqs. (18)-(21), we will have a 
potential V(x) — fA{x)(Jx, that can be eliminated from 
the dynamical equation for ip through the transformation 
-0 = e~*^''('^^(/) that turns Eq. (Ilj) into icf) = —iax(j>' + 
az mcj). Thus the potential can be exactly eliminated. 

Some deviations to the exact equivalences between free 
and potential-dependent dynamics presented here can be 
considered. For example, in the Dirac case a potential 
with/i(a;) = /3(x) = 0. 



In the cases in which F2{x) ~ 0, we can cancel the expo- 



V{x) = f2{x) (Jz + j/^{x) ax, 



(23) 



can be erased from the dynamical equation of ^ through 
the transformation 



where 



F^{x) = -if2{x), Fiix) = h{x), 



(24) 



(25) 



This mapping turns Eq. ([I]) into 



g ~iF2{x)<yy-iFi{x)^\ _ ^iF2{x)ay-iF4.(x) 

{-i(JxCJ)' +azm(l)). (26) 

Although we cannot eliminate the exponentials in both 



sides of the above equation as in previous cases, Eq. (26) 



is an approximate Dirac equation in regions where 
^2(0;) ~ 0. Notice that this does not necessarily entail 
/2(x) ~ 0, since the relation expressed in Eq. (251 shows 



that we have a condition between F2{x) and f2{x) 



D. Massless system 

We show now how to include a Dirac mass term into 
a massless dynamics using previous techniques. For 
massless particles the Majorana and Dirac equations are 
equivalent and can be written as 



i^jJ — —iaxTp + Vip. 



(27) 



We consider now that the potential in the above equa- 
tion has the form V{x) ~ /2(a^)cz- Through the trans- 
formation ■0 = e^^^^'^^^'^^cj) and imposing the constraint 
F!2{x) = -if2{x), Eq. (pTl) becomes 



nentials in both sides of Eq. ( 28 1, obtaining a free dynam- 
ics for the wave function (p. This condition is fulfilled by 
F2{x) = —imx when a; ~ 0, such that the transformation 
-0 — e^"^^"'y(f> relates the dynamics between massive and 
massless particles. That is, the dynamics of Eq. (27) with 
V = cr^TTi (/2 = m) is mapped onto the free dynamics in 
a small region near the coordinate origin. 



E. Examples 

We illustrate now our method with some examples. 
First, we will consider the ID Dirac and Majorana equa- 
tions with a linear potential. 



V{x) ^ gx, 
and the transformation, 

ip[x,t) — e 2 



ix,t). 



(29) 



(30) 



In this case, <j) satisfies the free Majorana equation and 
I0P — I0P, such that a solution of the Majorana equa- 
tion with potential has the same probability density of a 
solution of the free Majorana equation. The same trans- 
formation does not work in the Dirac case, shedding light 
on the different behavior of Majorana and Dirac particles 
under such potentials jll) . 

In Fig. [1] we show that our method works as a good 
approximation in a certain region of space in the case of 
the Dirac equation with a potential. 



V{x) — g cos(Ax) a^, 
and the corresponding transformation, 



Tp{x,t) 



'(j){x,t). 



(31) 



(32) 



-iF2{x)(7y 



JF2(x)ay 



lUx 



(28) 



Notice that the free dynamics (Fig. [TJd) of the trans- 
formed state is very similar to the dynamics under the 
potential of the untransformed state (Fig. flk) even in 
regions of space where the effect of the potential is non- 
trivial -as can be seen by comparing Fig. Ilk with Fig. 
[1];, where the free evolution of the untransformed state 
is plotted. 



III. BIPARTITE SYSTEMS. 

Our previous results can be extended to two-particle 
systems. For instance, it has been shown that a Lorentz- 
invariant two-body Dirac equation with an oscillator-like 
interaction can be written in the center of mass reference 
frame |15) . In ID and our particular representation, we 
have 

iip = 1= {ai—a2){ip'+mujxl3i2i>) + {l3i+l32)m'ip, (33) 

v2 
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FIG. 1. (Color online) a) Evolution of the probability density of a Dirac fermionic wave packet under the potential of Eq. 
(31 1, b) evolution of the corresponding state with the transformation in Eq. (321 under a free Dirac equation, c) evolution for 



the free Dirac case without transformation. For this example, we have chosen m = A, g — 2, and A — 15. 



with ai = (T^ (g) 1, a2 — 1 <E) a-^, Pi — a^ ® 1, h — 
1 ® (Tz, /3i2 = <yy® (Jy, and x — {xi — X2)/\/2- Thus, the 
corresponding two-body Majorana oscillator equation is. 



iV" = 7= («! - ^2)0^' + mu}x(3i2ip) - i{Pi + (32)m^*, 

v2 



with /3i = (Ty (g) 1 and $2 



explained above, we find that the mapping 



(34) 
ay. With the techniques 



tp{x,t) 



{x,t), 



(35) 



transforms Eq. ( |34| into a free two-body Majorana equa- 
tion near the coordinate origin. 



■|3l2 



■/3l2 



V2 



{ai~a2)(j}'~ii$i+f32)m cj)* 

(36) 

Thus, in this case the situation is similar to the one- 
particle example of Fig. fl] i.e., Eq. (34) can be simulated 
with good approximation with a two-body free Majorana 
equation in a certain region of space. Additionally, the 
probability density is not conserved under the transfor- 
mation. Interestingly, the same does not work for the 
two-body Dirac oscillator. 



the laser Rabi frequency, a^^-^x,z are the corresponding 
spin operators, and a, a^ are the phononic annihilation 
and creation operators. For general transformations, i.e., 
arbitrary functions of x coordinate times Pauli matri- 
ces, one may generate a unitary operator U that imple- 
ments an arbitrary continuous-variable state by concate- 
nated application of Jaynes-Cummings plus carrier inter- 
actions ^7\. To obtain the coupling to the spin degree 
of freedom, standard quantum conditional logic tech- 
niques for implementing the corresponding controlled-C/ 
gate may be used |18j . In general, with our formalism, 
the local transformation can be non-unitary, but always 
will map a pure state onto a pure state, given that it is 
just a matrix operator acting upon a vector state. This 
can be always realized by normalizing the final state and 
checking for the appropriate unitary transformation that 
connects initial and normalized final states. A normal- 
ization constant does not introduce any modification to 
the protocol given that all the considered equations are 
linear. In order to implement this kind of local transfor- 
mation upon a specific implementation like trapped ions, 
coupling between spin and motional degrees of freedom 
is provided by red and blue sideband laser pulses. A wide 
variety of static transformations can be implemented in 
this way [16j . 



IV. IMPLEMENTATION 

In order to realize this kind of protocol, the most in- 
volved part is the creation of the initial free-evolving 
state by the static transformation that encodes the po- 
tential dynamics onto free-evolving states. This can be 
much easier than the implementation of the correspond- 
ing potential in many different cases [161 117) . For exam- 
ple, in the case of the transformation in Eq. (301, this 



can be implemented in trapped ions by a straightfor- 
ward Jaynes-Cummings plus ant i- Jaynes-Cummings in- 
teraction in the dispersive limit, through the Hamiltonian 
H ^ hrjflia+e'^* + a-e-'^*){a + a'') ~ h[{T]nf/S]a,{a + 
a^)^ ex. x'^a^, plus a local rotation to change from the cr^ 
to the (Jx basis. Here rj is the Lamb-Dicke parameter, fl is 



V. CONCLUSIONS 

We have shown that Majorana and Dirac dynamics 
with a certain class of potentials can be simulated with 
the corresponding free dynamics by means of a static 
mapping between free states and states under the action 
of a potential. The Dirac and Majorana potentials for 
which this procedure is valid are of different nature, il- 
lustrating the different behavior of Majorana and Dirac 
dynamics [TT]. We also obtain that the probability den- 
sities of free and non-free states coincide for arbitrary 
scalar potentials V{x) = fi{x) in the Majorana equation. 
This implies that Majorana particles behave under this 
kind of potential as if they were free. In other cases, the 



method works as a good approximation in certain regions 
of space or parameter ranges. In particular, for massless 
Dirac systems a mass-like potential can be simulated in 
this way. We have extended our results to two-particle 
interacting systems. 
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